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Abstract 


The gravitational waves from the binary neutron star (BNS) merger GW170817 were accompanied by a multi-wavelength electromagnetic 
counterpart, which confirms the association of the merger with a short gamma-ray burst (SGRB). Afterglow observations implied that the 
event was accompanied by a narrow, ~ 5°, and powerful,10°° erg, jet. We study the propagation of a Poynting flux-dominated jet within the 
merger ejecta (kinematic, neutrino-driven and MRI turbulence-driven) of a neutrino-radiation-GR-MHD simulation of two coalescing neutron 
stars. We find that the presence of a post-merger low-density/low-pressure polar cavity, that arose due to angular momentum conservation, 
is crucial to let the jet break out. At the same time the ejecta collimates the jet to a narrow opening angle. The collimated jet has a narrow 
opening angle of ~ 4—7° and an energy of 10*? — 10° erg, in line with the observations of GW170817 and other sGRBs. 


Simulation 
“* The simulation for the BNS merger (Fig. 1) is — ait 
performed using a code described in Kiuchi et al. 
(2022); 106 
“* The NS is modelled with the SFHo equation of state 
(Steiner et al. 2013); 10° 
“ The binary is composed of 1.2 M, and 1.5 M, NSs. 
“* The corresponding chirp mass is consistent with the 10° 
one observed in GW 170817. 
“* The simulation domain is composed of 13 levels of ° 
the fixed mesh refinement (FMR) with the finest grid - 
spacing of Ax = 150 m and size of 
L € [—37.875,37.875] km; ; 
1 1 
“+ The merger and the subsequent black hole formation temp as 
take place at ~ 0.015 s and % 0.032 s, respectively. Fig 1. Left: Density map in a logarithmic scale of a 2D slice (x-z plane) through the center of the simulation box at t= 1s. We 
The formed BH has a spin parameter ofa = 0.65; superposed the velocity streamlines (white colour). Right: The jet parameter L calculated for the 2D slice of the ambient density 
¢* Total duration of the simulation is ~ 1.06 s for a jet luminosity of 10°” erg/s. 


Theory & Analysis 


“* As the jel propagates through the expanding ejecta it dissipates its | L, ) 
energy in shocks that form at the jet head. Jets expand time = 1.0 s » Mersgis] 
uncollimated if L = L,/X,p,c° > 1 (Fig. 1); 


4000 - 
“* If the jet is not powerful enough it may not break out from the 


ejecta (choked); 
¢¢ Escape condition E.., > E. 


jet yet, esc 


(from Gottlieb & Nakar 2022): 2000. 
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“¢ Collimation condition (Lyubarsky 2009): 
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¢¢ We calculate the contours with Eq. 2, assuming y = R/R, , we fii 

the contours and we find the relative local slope @ along the fit at ae 47.5 

the maximum height z, of the polar funnel, where the jet is free to “on =AnGn ab : a apa EHO 

expand (L > 1). (Fig. 2) x [km] 
«* Using the collimation contour and varying the luminosity of the jet 1020 1022 1024 1026 1028 1030 0 2 4 6 g 10 12 

L,, we are able to determine the relation between the luminosity P [erg cm~?] slope [deg] 

and a desired jet opening angle 0. We calculate the collimation Fig 2. Left: Map of the thermal pressure superposed onto the collimation contours that satisfy Eq. 2. The bottom color bar shows the values of the 

energy as E. Se ambient pressure while the different colours of the contour lines represent different jet luminosities L; (vertical color bar in log scale). We fixed the 

J,coll Ue light cylinder radius at R,; = 40 km. Right: Same as the left panel but superposed onto a power law fit z = ax’. The horizontal color bar at the 

«* We assume no delay time between the BH formation and the jet bottom of the figure represents the local slope of the curve fit (in degrees). 

launch in Eq. 1 (t7 = O) and engine times of 1, = 0.1 and 1 s (Fig. 3); t=0.3335 t=0.6675 t=1.000s5 
x i : 
x —— pay that manage to break out from the ejecta need Ea. t= barots Eapisaie 

J,eSC j + “~j,coll ° Fi cot, te =O.1 5 Fi col, te =O.15 Fi col, te =O.15 


Conclusions 


Fj, esc, magnetized F;, esc, magnetized F;, esc, magnetized 

“ The values we obtain (E, = 10°->*, 6, ~ 4° — 7°) are 
consistent with the observation of the afterglow light 
curve in combination with the superluminal motion of 
the jet core which estimate a powerful (E, = 10°-°° 


erg), and narrow (0, < 5°) conical jet; 


“* the jet can escape only at the poles; 


“* the engine operation time cannot be too short as 
compared to the time delay between the onset of the 
jet and the merger; ie eee nes 


Re Post-BN S merger ejecta anisotropic structure creates Fig 3. Escape energy, E, ..., (blue lines) and collimation energy, E, .., (red and green lines) versus 0 for three different snapshots at: 

t = 0.33 s (left), t = 0.667 s (center), and t=1 s (right). The red (and green lines) are drawn varying the luminosity of the collimation 

a favourable environment to successfully collimate contour and calculating, from the contour fit, the corresponding opening angle 0 measured at z, = 1670 km (left), z, = 3300 km (centre), 

jets in the energy range seen in SGRBs. and z, = 5000 km (right), which approximately represent the maximal heights of the low-density, low-pressure cavity at each snapshot, 
respectively. Solid blue lines represent the escape energy for magnetised jets. 


